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Abstract 

In various supersymmetric extensions of the Standard Model there appear non- 
topological solitons due to the existence of U{1) global symmetries associated with 
Baryon and/or Lepton quantum numbers. Trilinear couplings (A-terms) in the 
scalar potential break explicitly the U{1) invariance. We investigate numerically 
the stability of these objects in the case that this breaking is small. We find that 
stable configurations, oscillating though, can still appear. Other relevant properties 
are also examined. 
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Non-topological solitons are non-dissipative solutions of finite energy which arise 
in field theories possessing continuous global symmetries Of particular 

interest are possible stable configurations which carry baryonic or leptonic charge 
and appear in supersymmetric extensions of the Standard Model. In this latter 
context, they appear to have interesting cosmological consequences since they could 
be related to the baryon number of the universe, or even they could play a role as 
dark matter candidates 0] . In supersymmetric theories, the trilinear superpotential 
couplings of the matter and Higgs superfields, as well as the supersymmetry breaking 
terms, generate a scalar potential which may possess exact or approximate global 
symmetries which could be correlated to the above quantum numbers. It is the 
purpose of this letter to examine a generic form of such a scalar potential and 
investigate the properties of the solitonic configurations. 

In general, the resulting scalar potential in a supersymmetric theory depends 
on various scalar fields, U = U{^i,^2, ■ ■ ■ ,^n), however, in the present work we 
will concentrate on the simplest case of only one field $, i.e. U = f/($), while the 
generalization is straightforward. The contributions from supersymmetry breaking 
and non-renormalizable terms in the scalar potential have the generic form [| 

/ |<T)|2(n-l) (hn (f)*n \ 

In (|l]), rris is a soft mass term, of the order of the supersymmetry breaking scale. 
Xnr is a dimensionless constant associated to the Yukawa coupling of the corre- 
sponding non-renormalizable term of the superpotential. A is the coefficient of the 
trilinear term in the scalar potential (the so-called A-term), while M is of the order 
of the unification scale. It is easily observed that in the case of a scalar potential 
with vanishing A-terms the global U{1) symmetry of the Lagragian survives. Then, 
for appropriate values of the soft supersymmetry breaking parameters and Yukawas 
which determine the coefficients of the various potential terms, a stable solitonic 
configuration does appear. However, the A-terms are usually non-vanishing. Even 
if their initial value is zero, renormalization group running will generate a non- 
vanishing value and the aforementioned global U{1) symmetry breaks. In this letter 
we wish to investigate the stability of the solitonic configuration in the presence of 
small non-zero A-terms which induce small perturbative terms in the scalar poten- 
tial. 

To this end, we will use a simplified model in 1 -|- 1 dimensions and adopt the 
conventions and notations of reference [^. Assuming a complex scalar filed $, the 
dynamics of the resulting scalar field theory can be described by a Lagrangian density 

^ see for example @ , . 
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of the form 

C = ^df,<l>d^<^* - f/($) (2) 
Starting with a potential respecting a U{1) invariance 

f/($) = im2|$|2-ia|$|3 + i5|$|4 (3) 

while making use of the rescalings ||^ 

$ — $, X — (4) 
a m 

one obtains the Lagrangian 

2 ^ 2' ' 3' ' 4 ' ' ^ ^ 

with the parameter B = . Note that in order to get the above Lagrangian we have 
divided by a factor /a? which has dimensions of m? . Therefore, the Lagrangian 
in is dimensionless. The equation of motion reads 

+ |<|,|$ + 5|$|2$ = o (6) 

Using the ansatz [|| 

^{x) = a(a;)e'"* (7) 
and inserting (J^) in (P) the following differential equation for a{x) is obtained 

a" + {uj^ - l)a + a^-Ba^ = (8) 

The requirement for a finite energy configuration and the asymptotic behaviour at 
infinite "time", x oo, imply the following conditions for a: 

a'{0) = 0, ct(oo) = 0. (9) 

The constraints imposed on u by these conditions are written as 

< 1 - ^ < < 1 (10) 

The conserved charge is 

Q = yJ (**'9o$ - ^do^*)dx (11) 
2 



Given the ansatz (|^ one finds that the conserved charge is given by 

Q=cu J a{xfdx (12) 



Finally, the energy density is 



which, with the use of (0), becomes 

8 = + 1(1 + oo^y - -a^ + -Ba' (14) 
2 2^ ^ 3 4 ^ ' 

A numerical analysis of the model shows that static field configurations are obtained 
for a large portion of the parameter space a;^, B. Finally, the extension in 2 + 1 
dimensions shows that similar static stable 'objects' do exist, at least for certain 
cj^-values g]. 

We wish now to add a U{1) breaking term in the (rescaled) Lagrangian (^ of 
the form 

+ (15) 

n'. 

and examine in a similar manner the stability of the above configurations. We choose 
to work with n = 2 and find that the equations of motion for the scalar field in this 
case generalize as follows 

+ $ -|$|$ + B|<l>|2 + 2s$* = (16) 
<|,*_<|,*" + $*_|$|$* + 5|<|,|2 + 2s$ = (17) 

The current is given by ^*d^^ — $9^$* and the charge Q cx /($*(9o$ — ^do^*)dx, 
respectively, however, due to the extra f/(l)-breaking term the charge is no longer 
conserved. Its rate of change is given by 



dt 



J ZS($2 - <t>*^)dx 



An analytic solution of the above system is hard to find, so our intention is to 
solve numerically the complete equations of motion, taking as starting point the 
'unperturbed' solution (^. 

We first check that our numerical methods reproduce the original results in 
the case of vanishing s-term ( the [/(l)-invariant Lagrangian). Indeed, we have 
checked that the soliton-profile (the absolute value of the solution), to a very good 
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Figure 1: The time evolution of the sohton profile, the oscillation of the points 
X = 0, 2, 5, 7 and the parametric plot (Real Part, Imaginary Part) for one period 
and for three values of the breaking parameter s = 0.007, 0.003, 0.001. 

approximation, remains constant in time while the charge, energy and energy/charge 
are constant within 2%, 1.4% and 0.02% correspondingly. 

Next, we switch on the symmetry breaking term in the potential. In our numer- 
ical investigations we use the values B = 0.34 and tu^ = 0.51 for the two potential 
parameters. In Fig.(|l|) we show the evolution of the profile in time, the oscillation of 
\^{x, t)\ in time for four points {x = 0, 2, 5, 7) (x = is the top of the profile), and 
finally the parametric plot (Real Part, Imaginary Part) for the top of the profile. 
The last one is drawn for one period of the beat-like graph shown in the middle 
figure. The three rows correspond to three values of s = 0.007, 0.003, 0.001. 

Some comments are in order. Although the U{1) is explicitly broken, we see 
that we have a stable configuration though oscillating in time. We have checked the 
stability for a time corresponding to at least 10 periods {2it/uj). All points of the 
profile oscillate in phase, while the "movement" of a point resembles a beat. The 
higher we are on the solitonic profile the larger the oscillation is. The parametric 
plot, which for a soliton is a circle, transforms to a "turning ellipse" which closes as 
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Figure 2: The real (left) and imaginary part of the solution as a function of time, 
for X = and s = 0.007. 

a beat period is completed. As s gets smaller, the oscillation decreases, the top of 
the profile becomes smooth with the time while in the parametric plot, the "turning 
ellipse" becomes a circle. The period of the beat remains constant (within numerical 
errors). The non-smooth ending of the profile for large \x\ is also due to numerical 
errors in the solution of the equation of motion. 

It is interesting to analyze the real and the imaginary part of the solution, for 
X = 0, as a function of time; this is shown in Fig.(|^). We see a clear oscillation with 
frequency u which is modulated by a lower frequency. We clearly see the cosine 
and the sine function corresponding to the real and imaginary part (one should 
bear in mind that the ansatz for the solitonic configuration is a{x) exp{iujt)). The 
modulating signal has half a period phase difference between real and imaginary 
part. It is easy now to see that the high frequency of the beat-like profile is 2uj. 

As the breaking parameter s gets bigger, the oscillations show large amplitudes 
and the sense of a period fails to appear. This loss of clear periodicity starts to 
appear when the ratio energy/ charge, which also oscillates with time, becomes larger 
than one, which is the constraint for a stable configuration against decaying to the 
fundamental mesons of the theory. 

At this point, we should explain what we mean by charge. In the U{1) invariant 
theory the conserved charge is defined by the Eq.([TT|). We continue to define the 
charge by the same equation, in the sense that our symmetry breaking parameter 
is small. The energy is given by Eq. (|TT|) , where of course we have added the term 
coming from the breaking, namely (l/6)si?e[$^]. The lack of periodicity is most 
clearly seen in the parametric plot (Real Part, Imaginary Part) where now the line 
tends to cover all the plane, exhibiting two "fixed points". In Fig(^ we show the 
parametric plot for x = 0, s = 0.015 and for t = 300. The "fixed points" could be 
attributed to a left over Z2 symmetry. Indeed, the symmetry breaking term $"-|-$*" 
shows a Zn symmetry if $ = crlx) exp{iujt). 

We have made a Fourier analysis of the real and imaginary part of the complete 
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Figure 3: Parametric plot (Real Part, Imaginary Part), for x = 0, s = 0.015 and for 
t = 300 



solution for X = and, as it is obvious from the Figures, we get large amplitudes just 
for two frequencies which are multiplets of the basic frequency {2tt/T, where T is the 
period of the beat) namely u, which has the largest amplitude and u' = uj±(27i/T). 
To a very good approximation we can write 

$(x = 0, t) = cr cos(ti;t) + c'^ cos^uj't) + i{ci sin(co't) + c'j sin(u;'t)) (19) 

where the c's are constant coefficients. Since all points oscillate in phase we can 
write for the complete solution 

$(x,t) = K(x)$(0,t) (20) 

where K{x) is an unknown function of the position. The first obvious choice for 
K{x) is cr(x) itself. We have tried therefore the solution 

^x,t) = ^mt) (21) 

(T(0) 

which for the central region of the profile differs from the numerical solution by as 
low as 2%. Fourier analysis can also help to understand the transition to the simple 
solution, Eq.(|^). As s ^ 0, the amplitude of the u' frequency gets smaller and the 
only surviving frequency is the u. 

The above analysis holds when the parameters B and u of the solution are 
deep inside the stability region of the (i?,a;^) space. The stability region is defined 
from the Eq. ([lO|) [If. When we are near the boundary, the situation gets extremely 
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complicated. Even with a relatively small breaking parameter s, the energy/charge 
ratio becomes larger than 1. The period of the beat is no more constant but depends 
on the value of s. The Fourier analysis of the solution shows that more than two 
frequencies have large amplitudes and, for large enough s, none of them is equal to 
u. As s — 0, all the amplitudes, except one, tend to zero, while the remaining one 
corresponds to a frequency which tends to ui (since the beat period is not constant, 
the basic frequency of the Fourier analysis, and therefore its multiplets also change). 

We get a more clear situation near the boundary of the stable solutions, in the 
case of the so called thin wall approximation: the function a{x) has a constant 
value for a certain region in x and zero elsewhere. In that case, the Euler-Lagrange 
equation, the current and the energy density of the soliton are given respectively: 

$ + $ _ |<|,|$ + 5|<|,|2$ + 2s$* = (22) 
f = ^($*$ - $$*) (23) 

£ = l|<i.|2 + i|$|2 _ i|$|3 ^ :5|$|4 ^ S 2 ^ ^*2N ^24) 

2' ' 2' ' 3' ' 4 ' ' 2^ ^ ^ ^ 

The absence of spatial derivative terms in the energy makes the soliton more stable. 
We have found that on the boundary, the beat frequency is proportional to the only 
free parameter of the (rescaled) Lagrangian, namely B, taking care always to keep 
the breaking parameter s low enough to avoid the energy/charge ratio to become 
greater than 1. We have also checked that the same behaviour persists when we add 
a $^ term in the U{1) preserving potential. 

A last point to mention is the tremendously different situation we are facing 
in the case we change the power of the symmetry breaking term Eq.(|15[) in the 
Lagrangian to n = 3. In Fig.(^) we show the parametric plot for the same values of 
the parameters, B = 0.34, u"^ = 0.51 and for s = 0.05. Although the graph seems 
chaotic with respect to the corresponding one in the previous case, we clearly see 
three "fixed points". A larger s value is also needed, with respect to the previous 
case again, in order to see a significant disturbance of the soliton profile. There 
appears again a $(x, t) which is periodic in time, but it does not show a beat-like 
shape and the period seems smaller than before. We are planning to make a more 
general presentation of these points in a forthcoming publication. 

In conclusion, motivated by the U{1) breaking symmetry terms that appear in 
the soft-term Lagrangian, we investigated the simple case of a complex scalar field, 
in 1+1 dimensions, possessing a non-topological solition-type solution, with an extra 
explicitly breaking term. We solve numerically the equation of motion and found 
an oscillating with time qball-like solution. The points of the profile oscillate with a 
beat-like movement showing stability in time. As the breaking parameter gets larger 
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Figure 4: Parametric plot (Real Part, Imaginary Part), in the case where the break- 
ing term is $^ + and for x = 0, s = 0.05 and for t = 100 



and the oscillating energy/charge ratio gets greater than 1, instability appears in 
the solution in the sense of non constant period and with the destruction of the 
beat-like shape. 



We would like to thank P. Dimopoulos, K. Farakos, A. Kehagias and G. Kout- 
soumbas, G. Tiktopoulos for helpful discussions. 
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